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INTRODUCTION 


This  report  summarizes  activities  related  to  ONR  grant  N00014-97-1-0879  to 
the  Georgia  Tech  Research  Corporation,  for  the  period  1  June  1997  through 
September  30  1997  (Phase  I).  Ilie  work  during  Phase  I  has  focused  on  (1)  a 
sensitivity  analysis  of  the  effects  of  irrigation  on  the  biogeochemistry  of  coastal 
marine  sediments,  2)  the  upgrading  of  the  existing  multicomponent  reactive 
transport  model  for  early  diagenesis,  and  (3)  the  preliminary  development  of  a 
stochastic  algorithm  to  represent  burrowing  activities  of  benthic  macrofaima. 


ACHIEVEMENTS 


Reactive  transport  code 

The  basis  for  the  modeling  efforts  is  STEADYSEDl,  a  multi-component 
reactive  transport  code  for  steady-state  early  diagenesis  which  fully  incorporates  the 
reaction  couplings  among  the  reactive  species  of  C,  O,  N,  S,  Fe  and  Mn.  The  model 
has  been  tested  against  extensive  data  sets  collected  by  Canfield  et  at  (1993a,  b)  at 
three  coastal  marine  sites  that  exhibit  high  rates  of  combined  iron  and  manganese 
(hydr)oxide  reduction.  Results  show  that  the  model  provides  a  consistent 
explanation  of  the  entire  body  of  multi-component,  multi-site  observations.  The 
measured  concentration  profiles  of  28  individual  pore  water  and  solid  sediment 
species  are  satisfactorily  reproduced.  Furthermore,  the  model  predicts  the  observed 
distributions  of  sulfate  reduction  rates,  as  well  as  diagnostic  features  of  the  pore 
water  pH  profiles. 

The  coupled  nature  of  the  reactive  species  distributions  constrains  the  values 
of  reaction  and  transport  parameters,  because  a  single  set  of  parameters  must 
account  for  the  distributions  of  all  the  species  at  a  given  site.  The  parameters  are 
further  separated  into  site-dependent  (e.g.,  deposition  fluxes,  bottom  water 
composition,  mixing  coefficient,  porosity)  and  reaction-specific  parameters  (e.g., 
rate  coefficients,  apparent  equilibrium  constants,  limiting  substrate  concentrations). 
When  fitting  model-calculated  distributions  to  observed  ones,  care  is  taken  to 
minimize  variations  of  the  reaction-specific  parameters  from  one  site  to  another. 
This  modeling  strategy  emphasizes  extraction  of  reaction  parameters  from  field 
data  that  are  mechanistically  meaningful. 

Effects  of  irrigation 

A  sensitivity  analysis  was  performed  with  STEADYSEDl,  under  boundary 
conditions  representative  of  coastal  marine  depositional  environments. 
Conditions  were  chosen  on  the  basis  of  previous  work  (Van  Cappellen  et  al.,  1993; 
Tromp  et  al.,  1995;  Middelburg  et  al.,  1997).  The  results  show  that  model-calculated 
benthic  exchange  fluxes  of  most  solutes  are  dominated  by  irrigation.  Because  of 
high  rates  of  molecular  oxygen  reduction  in  organic-rich  coastal  sediments,  the 
contribution  of  molecular  diffusion  to  benthic  uptake  of  O2  may  in  some  cases  be  of 
the  same  order  of  magnitude  or  even  larger  than  that  of  irrigation.  For  nitrate. 


1 


molecular  diffusion  and  irrigation  cause  exchange  fluxes  in  opposite  directions.  As 
a  consequence,  there  is  a  net  transfer  of  nitrate  from  the  water  column  to  the 
sediments,  although  the  interfacial  pore  water  gradients  predict  diffusional  fluxes 
out  of  the  sediments.  The  results  highlight  the  fundamental  importance  of  pore 
water  irrigation  by  macrofauna  in  controlling  solute  exchanges  between  sediments 
and  water  column  in  coastal  environments. 

Benthic  fluxes  and  chemical  conditions  in  sediments  not  only  depend  on  the 
integrated  intensity  of  irrigation,  but  also  on  the  depth  distribution  of  the  irrigation 
exchange  coefficient.  This  is  shown  in  Figure  1,  where  three  different  distributions 
of  the  vertical  irrigation  exchange  coefficient,  a  (in  units  of  yr^),  are  considered.  In 
two  cases  (1  &  3),  the  irrigation  coefficient  is  assumed  to  decay  exponentially  with 
depth.  This  is  the  distribution  used,  for  example,  by  Martin  and  Banta  (1992)  to 
reproduce  tracer  distributions  of  pore  water  bromide  and  222^^/226^3 
disequilibrium  in  Buzzards  Bay,  Massachusetts.  In  the  other  case  (2),  the  depth 
distribution  of  a  follows  a  step  fimction,  that  is,  a  uniformly  irrigated  surface  layer 
overlies  non-irrigated  sediment.  The  depth  integral  of  a  is  the  same  for  all  three 
distributions.  Thus,  the  simulations  compare  sediments  experiencing  the  same 
overall  irrigation  intensity. 

The  calculated  pH  profiles  show  distinct  differences  between  the  three 
irrigation  scenarios  (Figure  1).  (Note:  pH  is  one  of  the  most  sensitive  parameters  of 
sediment  biogeochemistry.)  The  differences  are  most  pronounced  when 
comparing  the  step  function  to  the  exponential  distributions.  In  the  depth  range  2- 
3  cm,  the  predicted  pH  values  may  deviate  from  one  another  by  as  much  as  one  pH 
unit,  i.  e.  the  equivalent  of  an  order  of  magnitude  in  hydrogen  ion  concentration. 

Figure  1.  Irrigation  scenarios:  (a)  imposed  depth  distributions  of  the  irrigation 
exchange  coefficient;  (b)  calculated  pH  distributions.  Calculations  were  performed 
vmder  conditions  typical  of  a  coastal  marine  sediment.  Boundary  conditions  and 
parameter  values  are  the  same  in  the  three  simulations,  except  for  the  distribution  of 
the  irrigation  exchange  coefficient.  See  text  for  further  discussion. 


Irrigation  coafficiant  (1/a)  pH 

^  0  100  200  300  400  6.5  7  7.6  8  8.6  9 

0 


2 


Existing  early  diagenetic  models  use  idealized  depth  distributions  of  the 
irrigation  coefficient  such  as  those  presented  in  Figure  1.  Nonetheless,  little  data  or 
theoretical  arguments  have  been  presented  for  selecting  one  distribution  over 
another  (for  exceptions,  see:  Martin  and  Banta,  1992;  Wang  and  Van  Cappellen, 
1996).  The  sensitivity  analysis  performed  here,  however,  shows  that  the  depth 
distribution  of  the  exchange  coefficient  has  a  significant  impact  on  sediment 
biogeochemistry.  A  large  part  of  the  work  planned  for  Phase  n  of  the  project  will  be 
devoted  to  the  development  of  better  theoretical  and  experimental  constraints  on 
the  functional  dependence  of  a  on  depth  below  the  water-sediment  interface. 

CaCOs  and  pH  buffering 

Dissolution  of  deposited  biogenic  CaCOa  or  precipitation  of  authigenic  CaCOs 
may  exert  a  critical  control  on  pore  water  pH,  alkalinity  and  total  dissolved 
inorganic  carbon  (e.g.,  Gaillard  et  al.,  1989;  Cai  et  al,  1995).  The  published  version 
of  STEADYSEDl  (version  1.1)  does  not  include  the  dissolution  and  precipitation 
reactions  of  CaCOs.  In  the  absence  of  buffering  by  these  reactions,  the  model 
simulations  tend  to  overpredict  the  magnitude  of  pore  water  pH  fluctuations. 

STEADYSEDl  was  upgraded  (version  1.2)  by  including  the  rates  of  CaCOs 
dissolution  and  precipitation  in  the  continuity  equations  of  Ca2+,  CaCOs,  total 
inorganic  carbon  (TC)  and  total  alkalinity  (Aik).  As  before,  the  pH  profile  is 
calculated  by  assuming  that  the  dissolved  carbonate  system  remains  in  local 
equilibrium.  This  is  justified  by  the  fast  homogeneous  acid-base  kinetics.  Several 
rate  models  for  CaCOs  available  for  the  user  to  choose  from,  as  summarized  in 
Table  1.  Preliminary  results  show  that  by  explicitly  including  CaCOs  dissolution, 
the  pore  waters  are  buffered  against  large  decreases  in  pH.  More  simulations  are 
needed  to  fully  assess  the  effects  of  CaCOs  precipitation. 


Table:  Options  for  CaCOs  dissolution/precipitation  kinetics  in  Steadysedl  (v.  1.2) 


Model 

Rate  expression 

Reference 

Empirical 

Rj  =  kj[CaCOs]{l-^lf^ 

Rd-.  Cai  et  al.  (1995) 

Rp.  Zhong  and  Mucci  (1993) 

Mechanistic  I 

+  k2au^co^  + 

Rp=k.,ac^,.a^co-, 

Chou  et  al.  (1989) 

Mechanistic  11 

Real  =  k^lH*]  +  k2lH^C0U-i-  k^lH^O] 
-k^lCa^^llHCO;] 

Plummer  et  al.  (1979) 

Mechanistic  III 
(surface  reaction) 

Ratal  =ki>  CO^al^  +  (fcj  -  fcj)  >  Ca^a„^co, 
+k,-(k,-k,)>CO,H\^^^^. 

-k-j  >  Ca*a„^co^a^aCol  ~  ^^caCO° 

Arakaki  and  Mucci  (1995) 

R<j:  dissolution  rate;  Rp.  precipitation  rate;  Rtotai-  overall  reaction  rate;  Ui:  activity  of 
solute  i;  [i]:  concentration  of  solute  f;  kf  rate  constant;  >i:  density  of  surface  species  i 
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The  new  version  of  the  reactive  transport  code  has  not  yet  been  made 
publicly  available  because  some  numerical  problems  remain  to  be  solved.  These 
problems  are  the  consequence  of  the  fast  rates  of  CaCOa  dissolution  and 
precipitation,  relative  to  the  other  kinetic  processes  represented  in  the  model  (e.g., 
the  degradation  of  organic  matter).  Large  differences  among  rate  terms  in  the 
continuity  equations  are  prone  to  cause  numerical  instabilities  in  the  finite 
difference  solutions.  We  are  currently  implementing  changes  in  the  code  that 
counter  these  instabilities. 

Stochastic  burrowing 

First  steps  have  been  taken  in  the  development  of  a  two-dimensional  (2-D) 
stochastic  model  for  random  vertical  burrowing.  In  this  model  the  water-sediment 
interface  is  represented  by  a  set  of  equally  spaced  spatial  nodes.  The  algorithm 
considers  each  node  in  succession  and  decides  whether  or  not  to  initiate  burrowing. 
The  decision  is  based  on  a  randomly  selected  number,  whereby  the  relative  range  of 
values  leading  to  burrowing  is  determined  by  the  average  normalized  burrow 
density. 

The  burrowing  strategy  plus  environmental  constraints  that  direct  the 
activity  of  the  macrofauna  are  translated  into  a  probability  distribution  for  the 
burrow  length.  Once  a  node  has  been  selected  as  the  site  of  a  burrow,  the  length  is 
determined  by  generating  another  random  number.  The  value  ranges  assigned  to 
different  lengths  are  proportional  to  the  corresponding  probabilities.  Based  on 
numerous  realizations  of  the  stochastic  burrow  model,  we  calculate  the  average 
depth  distribution  of  available  burrow  surface  area  (a  proxy  for  the  exchange 
coefficient),  plus  the  deviation  from  the  average.  In  this  manner  a  statistical,  rather 
than  a  deterministic,  description  of  sediment  burrowing  is  achieved.  Our  efforts  in 
the  near  future  will  focus  on  the  formulation  of  realistic  probability  distributions 
for  burrow  depths. 


OTHER  ACnVTTTES 

A  meeting  of  the  P.I.s  and  the  Program  Officer  was  held  at  the  Naval 
Research  Laboratory  in  Stermis  on  September  26,  1997.  Progress  to  date  by  the 
different  laboratories  was  reviewed,  and  plans  for  Phase  11  (12/1/97  to  9/30/99)  were 
discussed.  Emphasis  was  placed  on  strengthening  the  integration  of  the  P.I.s' 
efforts. 

Work  has  begun  on  a  manuscript  on  the  effects  of  bio-irrigation  on  solute 
exchanges  using  data  from  the  June  97  cruise  to  the  Dry  Tortugas  organized  by  NRL 
(Lavoie,  chief  scientist).  Dr.  Furukawa  is  the  lead  author  of  this  collaborative  effort. 

An  abstract  based  on  research  closely  related  to  that  supported  by  this  grant 
was  submitted  to  the  special  session  on  Diagenesis  in  Coastal  and  Marginal 
Sediments  of  the  1998  Ocean  Sciences  Meeting  (San  Diego,  February  9-13, 1998).  It  is 
entitled  "Spatial  and  Temporal  Variability  of  the  Iron  and  Sulfur  Cycles  in  Salt 
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Marsh  Sediments:  An  Integrated  Geochemical  and  Microbiological  Study",  by  P. 

Van  Cappellen,  A.  Roychoudhury,  E.  Viollier,  K.  Lowe,  T.  DiChristina,  T.  and  J. 

Kostka. 

The  following  activities  are  in  the  planning  stage;  they  will  enhance  the 

visibility  of  research  activities  supported  by  ONR  grant  N00014-97-1-0879. 

•  A  series  of  invited  lectures  by  the  P.  I.  at  the  Graduate  College  of  GEOMAR 
(Forschungszentrum  fiir  marine  Geowissenschaften)  in  Kiel,  Germany  (Nov. 
26-28,  1997).  The  lectures  will  focus  on  the  theoretical,  numerical  and 
experimental  aspects  of  reactive  transport  in  marine  sediments. 

•  The  organization  by  the  P.  I.  of  a  thematic  session  entitled  "Biogeochemical 
Dynamics:  From  Molecular  Approaches  to  Predictive  Models"  at  the  Conference 
on  Southeast  Coastal  Ocean  Research  (SECOR)  (Savannah,  April  7-10,  1998). 


PERSONNEL 


In  addition  to  the  P.  I.,  first-year  graduate  student  Chao  Liu  and  Postdoctoral 
Fellow  Eric  Viollier  have  been  involved  in  upgrading  the  STEADYSEDl  code  and 
running  sensitivity  simulations.  Meanwhile  Dr.  Viollier's  appointment  at  Georgia 
Tech  has  ended.  In  anticipation  of  continued  funding  for  Phase  II  of  the  project,  a 
new  Postdoctoral  fellow  has  been  recruited  to  work  on  the  stochastic  bio-irrigation 
model:  Dr.  Carla  Koretsky  (Ph.  D.  1997,  Johns  Flopkins  University)  will  join  the 
P.I.'s  research  group  on  January  1, 1998. 
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